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1 . Introduction . 

This publication has been compiled as the result of three 
courses given at the TOOL-ORP windmill project in (Jhazipur. 
Those courses lasted from 10 to 8 weeks and they were a 
combination of theoretical aspects and practical training 
in the design, construction and installation of the 
12 PU 5 CO windmill. 

Another part of the courses was the implementation of the 
windmills with irrigation aspects, economical aspects and 
social aspects. This "agricultural" part, however, is not 
included in this publication, but is published separately 
by TOOL/SWD and has been compiled by A.E.M. van Vilsteren 
(SWD 81-4). 

This publication aims at practical applications such as 
analyzing windregimes, determining rotor diameters, 
calculating bladeshapes, construction aspects of the 
windmill, pump choice, matching pump to windmill etc. etc. 
It is mainly aggravated to the 1 2 PU f>00 windmill. 
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2. WIND ENERGY 

2,1 Available windpower 

Wind is air in motion. If air with a mass m is moving with 
velocity v, it has kinetic energy expressed by: 

E = i mv 2 [j] (2.1) 

If the density of the air is p, then the kinetic energy 
per volume of air is: 

E v = *? y2 [ J /m 3 ] (2.2) 

For an area A perpendicular to the wind direction the flow 
per second through A is: 

V = v.A [m 3 /s] (2.3) 

So the power that flows with the air, through area A, is 
the product of the kinetic energy of the air and the air- 
flow: 

P = ifv 2 x v.A = £ov 3 A [w] (2.4) 

This is the power available in the wind. Only part of 
this power can actually be extracted by a windmill. 

The above derived relation for the power in the wind 
shows clearly that: 

- The power is proportional to the density p. 

This factor cannot be influenced and varies slightly with 
the hej/5ht and temperature. (For 20°C at sea level p e 
1.2 kg/m 3 .) 

- The power is proportional to the area swept by the rotor 
and thus to R 2 . Radius R is chosen in the design. 

- The power varies with the cube of the wind velocity v. 
Note that the power increases eightfold if the wind speed 
doubles. 
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Part of the kinetic energy of the airflow is carried over 
to the blades of a wind rotor resulting in mechanical power. 



3. AIRFOILS 

3.1 The relations between the velocity at the rotor blade and 
the forces acting on the blade. 

The rotor of a windmill consists of one or more blades attach- 
ed to a hub. The cross sections of these blades can have sev- 
eral forms; we call these cross sections air foils. An air foil 
is a surface over which air flows. This flow results in two 
forces: LIFT and DRAG. 

Roughly, lift is the force needed to bend the flow. 

Lift is always perpendicular to the airflow. 
Drag is measured parallel to the flow. 



UFT 




Lift and Drag 
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All air foils require some angle with the air flow in order 
to produce lift. The more lift required, the larger the 
angle (till a certain extent). 

The chord line connects the leading edge and the trailing 
edge of the airfoil. The angle required for lift is called 
angle of attack. The angle cf attack is measured between 
the chord line and the direction of the airflow. 




To describe the performance of an airfoil independent of 
size and airflow velocity, we divide lift L and drag D by 
^>v^A where 

p = air density |kg/m^] 
v = flow velocity [m/s] 

A = blade area ( = chordx blade length) [m 2 J 

The results of these divisions are called lift coefficient Ci 
and drag coefficient C, 




(2.5) 



(2.6) 



As stated before, the amount of lift and drag that is 
produced, depends on the angle of attack. 



This dependence is a given characteristic of an airfoil 
and is always presented in C-j-oH and C^- graphs. 

Cl 01/ 




Lift and drag characteristics 



The Cfl/C^ - ratio is choosen as small as possible in order 
to optimize the rotor design. 

For the design of a windmill the and ot values have to 
be found that correspond with a minimum - ratio. 

The following table gives these design values for some 
air foils. 
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Airfoil geometrical 
name description 


(c d / c i } 

min 


cx 


C l 


Sail and pole ^ ^ t 

£- — c J 


0.1 


5 


0.8 




0.1 


4 


0.4 


Arched steel plate J— r" " f/c=0.07 

\m c J f/c=0.1 


0.02 
0.02 


4 
3 


0.9 
1 .25 


Arched steel plate ^^^fl^^^ f/cs=0.07 
with tube on concave J L a 
side f/c=0.1 
Slde d<0.1 c 


0.05 
0.05 


5 
4 


0.9 
1 .1 


Arched steel plate ^^^f^^^^ 

with tube on convex f/c=0.1 

side 


0.2 


14 


1.25 



A blade element of a windmill rotor "sees" a relative 
velocity that results from the wind velocity in combination 
with the velocity with which the blade element moves itself. 




rotor plane 



tangential 
velocity 



Relative velocity on rotor blade. 

is the angle between the relative velocity W and the 
rotor plane. 
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4. WINDMILL CHARACTERISTICS 

4.1 Torque and power characteristics. 

The components in the plane of rotation of the lift forces 
result in a force working in tangential direction at some 
distance from the rotor center. This force is diminished 
by the component of the drag in tangential direction. The 
result of these two components is a propelling force in 
tangential direction at some distance from the rotor center. 

The product of this 
propelling force and its 
corresponding distance to 
the rotor center is the 
contribution of the blade 
element under consideration 
to the torque T of the 
rotor. The rotor rotates 
at an angular speed: 

-*1= 2TTn [rad/s] (4.1 ) 

The power that such a rotor 
extracts from the wind is 
transformed into mechanical 
power. This power is equal 
to the product of the torque 
and the angular speed. 

T s torque [nm] 
SI = angular speed [rad/s] 
P = T £L [w] (4.2) 

A windmill of given dimension transforms kinetic energy from 
the wind into a certain amount of power. Equation (4.2) clearly 
shows that a windmill for a high torque load (for example a 



L e Lift 




relative 
velocity W 



piston pump) will have a low angular speed; a high speed 
design will only produce a small amount of torque (for 
example for a centrifugal pump or an electricity generator). 

We call a graph that shows the dependance of the windmill 
torque on the angular speed, a windmill torque characteristic 
The following figures show the torque characteristics of two 
different windmills designed for the same power, taut for 
different angular speeds. The windmill torque characteristic 
depends on wind speed v so we have many curves in one charac- 
teristic. 



T 




windmill windmill 

With relation (4.?) it is very simple to derive from the 
torque characteristics the corresponding power chracteristics 
See following figure, where the power-angular speed curve, 
belonging to windmills of the above figure are shown. 
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low speed 
windmill 




high speed 
windmill 



Note 1 . The power of the two windmills is the same, but it 
is delivered at different angular speeds. 

2. The maximum power is delivered at a higher 
angular speed than the maximum torque. 

3. The maximum of the power curves vary with 
the cube of the angular speed: 

3 (4.3) 



max 



SI 



while the corresponding torque values vary 
with the square of the angular speed: 



T (at P = 



P ) 
max ' 



(4.4) 



The starting torque, i.e. the torque at zero 
revolutions per second, is considerably lower 
for high speed than for low speed windmills. 
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4.2 Power coefficient 

In order to be able to compare the properties and 
characteristics of different windmill design under different 
wind conditions we write the mechanical power as the power 
in the air multiplied by a factor C 



mech 



= C. 



air 



C p is called power coefficient and is a measure for the 
success we have in extracting power from the wind. With 
relation (2.4) we may write: 



C_ = 



mech 



h O V J TTR £ 



(4.5) 



For the same reason we divide the speed u of the rotor at 
radius r by the windspeed v. 




Definition of speed ratio 



We call the result local speed ratio which is noted 



X =r E 

r v 



SI. r 

v 



(4.6) 



Hie speed ratio of the element of the rotor blade at radius 
R is called tip-speed ratio: 



V 



(4.7) 



Note: later it will be shown that a windmill has one value 
of X at which the power coefficient is maximum. This 
is often called «the tip-speed ratio of a windmill » 
or »the speed ratio (Lambda) of a windmill* . 
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From relation (4.2) we know that 



(4.9) 



With this relation we define a dimensionless torque 
coefficient in the following way: 



P 



= C 5 f v 3 Tf R 



v 
"TT 



(4.5) 
(4.7) 



T = £ 



P 



(4.9) 



We define: 



1 

X 



If V^TT R j 



T 



c T = 



Ifv 2 TT R 3 



(4.10) 



Note that in this way relation (4.4) is still valid, but 
now in dimensionless form: 



(4.11 ) 



An empirical formula to derive the starting torque coef- 
ficient of a rotor as a function of its optimum tip speed 
ratio is: 

C _0.5 
T start ~ TT 
d 

^ d = design tipspeed ratio 
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4.3 Basic form of a windmill characteristic 

The power coefficient C p in equation (4.5) is not an 

efficiency, but may be interpreted as a measure of the 

success that a windmill has in transforming wind energy 

into mechanical energy. For one specific windmill C 

varies with tip-speed ratio of the windmill. In 

dimension-less form this is shown in a so-called C - X 

P 

characteristic based on formulas (4.5) and (4.7). 

See fig. below where one curve now represents all the 
curves for different v of a certain windmill. 




low speed high speed 




o i 2 3 *i S 6 J 1>\ 

Cp-A characteristic 



This characteristic is independent of air density 
windspeed v and radius R. 

Using relation (4.11) we may derive from the above fig. a 
dimensionless form of the torque-speed characteristic of 
the windmill: C T - X curve; see next fig. 



Also here one curve represents all curves of a certain 
windmill . 
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Ctp - X characteristic 



Note that the power is zero if X = 0, but that the torque 
is not. See relations (4.2), (4.11). 

4.4 Maximum power coefficient 

The power coefficient C as defined with relation (4.5) 

It 

describes how much power we get from the wind with a 
windmill. The power in the wind is given by relation 
(2.4). We are of course very interested in how much 
windpower we can transform into mechanical power with 
a windmill. In other words, we want to know what the 
highest power coefficient is for a given windmill that 
is designed for a certain tip-speed ratio. Betz was the 
first one to show that the theoretically maximum 
attainable power coefficient is 0.593- Three other effects 
cause a further reduction of the maximum power coefficient. 
How to find the maximum power coefficient that takes these 
effects into account will be explained in the following: 
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Betz coefficient 

It is not possible to transform all the wind energy that 
flows through cross sectional area A into mechanical energy. 
If we could transform all the energy in the air this would 
mean that we could extract all kinetic energy from the air; 
the air velocity behind the rotor would then be zero and no 
more air would flow through the rotor. The process of 
extracting kinetic energy from the wind will stop and no 
more power will be transformed. If on the other hand the air 
velocity behind the rotor is equal to the wind velocity, no 
kinetic energy has been extracted and also in this case no 
power will be transformed. In this way it may be understood 
that, if the flow velocity behind the rotor is either zero 
or equal to the wind velocity v, in both cases the mechanical 
power is zero. Between these values there is an optimum value of 
the wind velocity behind the rotor. 

Betz found this value to be -j v and calculated the 
maximum power coefficient (Betz coefficient). 



C 



P, 



= 0.593 



(4.12) 



max 



However, this value is only valid for a theoretical 
design for a high tip-speed ratio , with an infinite 
number of blades and a blade drag equal to zero . 
The effect of deviations from these three assumptions 
will be shown one by one in the next three paragraphs. 





'/ 3 v 



Maximum energy 
conversion for : 



A 
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Effect of wake rotation on maximum power coefficient . 

The Betz coefficient suggests that, independent of the 
design tip-speed ratio \ d , we may expect a maximum 
power coefficient Cp of 0.593- Relation (4.12) is, 
however, only valid for high tip-speed ratios and for low 
tip-speed ratios considerable deviations exist. This can 
be explained in the following way: 
The power is: torque x angular speed. The torque is 
produced by forces acting on the blades in tangential 
direction, multiplied by their corresponding distances 
to the rotor centre. These forces are the result of 
velocity changes of the air in tangential direction 
(action = reaction force = mass x velocity change per 
unit of time). The direction of the velocity change in 
the air is opposed to the direction of the forces acting 
on the blades. Since the air has no tangential velocity 
before passing the rotor, the velocity change means that 
behind the rotor the wake rotates in a direction opposite 
to that of the rotor. This wake rotation means a loss of 
energy because the rotating air contains kinetic energy 
(see relation 2.1). Since a certain amount of power is 
to be transformed, we know from relation (4.2) that a low 
tip-speed ratio (= low angular speed .a, means that the 
torque T must be high. High torque means large tangential 
velocities in the wake; the consequence is a loss of 
energy and a lower power; the more so if the design tip- 
speed ratio is lower. The result, shown in the graph 
presented below, shows the collection of maximum obtainable 
power coefficients of ideal windmills i.e. windmills with 
an infinite number of blades without drag. 



maximum power coefficients of ideal windmills 
Effect of C H /C-, -ratio on maximum power coefficient. 



The factor C^, as defined in relation (2.6) is a measure 
for the resistance of the blades against moving through 
the air. The Cd/Ci-ratio determines the losses due to this 
resistance. These losses are calculated and included in 
the collection of maximum power coefficients of the latest 
fig. The results are shown in fig. below. 




o I t 3 M 5 6 ? 8 9 to ii 12 IS \H IS — O X 

Effect of C^/C-j-ratio on Cp-raax for a rotor with 
an infinite number of blades. 
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This graph shows that a rotor, designed for a tip-speed 
ratio X = 2.5 with airfoils having, for example, a 
minimum Cft/Ci-ratio of 0.05, will have a maximum power 
coefficient C p = 0.46. If the rotor is designed for 
X = 10, at the same value, the Cp value atX = 10 

will have a maximum of 0.3. 

Note that from the same graph it is clear that it is 
useless to design a rotor for X= 10 with airfoils that 
have (C d /Ci) m in = 0.1 . 

Effect of number of blades on maximum power coefficient . 

The number of the blades also affects the maximum power 
coefficient. This is caused by the so-called "tip-losses" 
that occur at the tips of the blades. These losses depend 
on the number of blades and the tip-speed ratio. The los- 
ses have been calculated and as example are included in 
the collection of maximum power coefficients for (Cd/Pl^in 12 
0.03 of the above graph. The results are shown in graph below 




Influence of number of blades B on Cp-max 
for Cd/Cx = 0.03. 
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In the following graph the maximum expected power 
coefficient is given for X -design between 1 and 15, 
for a constant C^/C^ ratio of 0.02 (arched steel plate) 
while the number of blades B is varied. 

Conclusion: Tf design figures for tipspeed-ratio X and 
number of arched steel plate blades B have been chosen, 
the expected power coefficient may be read from this 
graph. 



Effect of number of blades for Cd/CI = 0.02 




5. DESIGN OF A WINDMILL ROTOR 

5.1 Calculation of blade chords and blade setting 

Earlier it was shown that the selection of the number of 
blades B affects the power coefficient. Although B has 
no influence on the tip-speed ratio of a certain wind- 
mill, for the lower design tip-speed ratios in general 
a higher number of blades is chosen (see following table). 
This is done because the influence of B on Cp is larger at 
lower tip-speed ratios. A second reason is that choice of 
a high number of blades B for a high design tip-speed ratio 
will lead to very small and thin blades which results in 
manufacturing problems and a negative influence on the lift 
and drag properties of the blades. 



Selection of 
number of blades. 



X 


B 






1 


6 




20 






2 


4 




12 






3 


3 




6 






4 


2 




4 


5 




8 


2 




3 


8 




15 


1 




2 



A second important factor that affects the power coefficient 
is the drag. Drag affects the expected power coefficient via 
the C d /C- t - ratio. This will influence the size and even 
more the optimal speed ratio (A ) of the design. 

A large C d /C 1 - ratio restricts the design tip-speed ratio. 
At lower tip-speed ratios the use of more blades compensates 
the power loss due to drag. See graph 1 . In this collection 
of maximum power coefficients it is seen that for a range of 



?1 

design speeds 1< A^1Q the maximum theoretically attainable 

power coefficients lie between 0.35 <iC <* 0.5 . 

P max 

Due to deviations, however, of the ideal geometry and hub 
losces for example, these maximums will lie between 0.3 
and 0.4. This result shows that the choice of the design 
tip-speed ratio hardly affects the power output. Two other 
factors, however, limit the choice of the design tip-speed 
ratio. One is the character of the load that in most cases 
will require a high starting torque, the design speed of the 
rotor will usually be chosen low; this allows the designer 
to use simple airfoils like sails or steelplates. If the 
load is fast running like a generator or a centrifugal pump, 
then a high design speed will be selected and airfoils with 
a low C^/C^ - ratio will be preferred. The second factor is 
that the locally available technologies will often restrict 
the possibilities of manufacturing blades with airfoils 
having low C^/C^ - ratios. But even in the case of high speed 
design, simple airfoils like arched steel plates, can give very 
good results. 

For design of the blades we need the following data: 



Rotor radius 


R [m] 


Humber of blades 


B [-] 


Design tip-speed ratio 


X [-] 


Airfoil data: Ca/Cj- ratio 


H 


Design lift coefficient 




Corresponding angle of attack 





Once these data are known it is now very simple to calcu- 
late the blade geometry; i.e. the chord c of the blade and 
the blade angle f3 , the angle between the chord and the plane 
of rotation. 
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Blade setting y9 



Only four simple formulas are needed: 



X 

r 


= X Q x r/R 


(5.1) 


c 


8TT r 

= O-co Sj 0) 


(5.2) 




= - <* 


(5.3) 





= | arctan ^ 


(5.4) 



The latter relation is also presented in table 1 , Appen- 
dix I. 

The underlying theory is too complicated to be explained 
here. The reader who is primarily interested in the design 
of the rotor, can do without this theory. 

Now the design procedure is as follows: 

Divide the blade with radius R in a number of parts of 

equal length. In this way we find cross sections of the blade. 

Each cross section has a distance r to the rotor centre and 

has a local speed ratio X r> according to (5.1). 

In appendix I we can find the corresponding angle for each 

cross section. is the angle of the relative air-velocity 

W that meets the blade section at radius r. We now calculate 

the chord with relation (5.2). For the sake of ease, (1-cos0) 
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has been added in the table of appendix I. The blade angle 
at, the corresponding radius is found with (5.3). 

As an example we will take the 12 PU 500 windmill with: 

R = 2.5 [m] 
B =12 

Airfoil: 10$ arched steel plate 
with: 

* 1.25 (i.e. the value for minimum (C^/C^) 
«* = 3° (i.e. the corresponding angle of attack) 

With equations (5.1) and (5.2), (5.3) we can now compute the 
values of the following table: 



Cross 

section 

number 


r[m] 


K 


0° 


*o 




C [m] 


1 


0.5000 


0.4 


45.6 


3 


42.6 


0.252 


2 


0.6250 


0.5 


42.3 


3 


39.3 


0.273 


3 


0.9550 


0.76 


35.2 


3 


32.2 


0.293 


4 


1.2500 


1.00 


30.0 


3 


27.0 


0.281 


5 


1.5625 


1.25 


25.8 


3 


22.8 


0.260 


6 


1.9100 


.1.53 


22.1 


3 


19.1 


0.235 


7 


2.1875 


1.75 


19.8 


3 


16.8 


0.217 


8 


2.5000 


2.00 


17.7 


3 


14.7 


0.198 



The result is the blade chord c and the blade setting p at 
various stations along the blades. Plotting the chords gives 
the blade form and plotting f& shows the desired twist of the 
blade. 
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cross section 12 34 5 67 8 



Blade chord and blade setting at various stations along 
the blade. 
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5.2 Deviation of the calculated chords and blade setting 

In the last paragraph we showed how to calculate the ideal 
blade form. The chords as well as the blade angles as cal- 
culated vary in a non-linear manner along the blade. Such 
blades are usually difficult to manufacture and lead to an 
uneconomic use of materials. In order to reduce these pro- 
blems it is possible to linearize the chords and the blade 
angles. This results in a small loss of power. If the line- 
arization is done in a sensible way the loss is only a few 
percent. 

In considering such linearizations it must be realized that 
about 75# of the power that is extracted by the rotor from 
the wind, is extracted by the outer half of the blades. 
This is because the blade- swept area varies with the square 
of the radius; also the efficiency of the blades is less at 
small radii, where the speed ratio X r is small. On the other 
hand, at the tip of the blade the efficiency is low due to 
the tip losses. 

For reasons mentioned above, it is advised to linearize the 
chords c and the blade angles fi between r=0.5 R and r=0.9 R 
without too much deviation to the computed theoretical values. 

Example: We linearize the blade chords c and the angles (5 as 
calculated . 

The nearest value of r to 0.5 R in the table is r= 1.25 
(= 0.5 R). The nearest value of r to 0.9 R is r= 2.1875 
(= 0.875 R). 

We can now linearize the chords by writing c in the following 
way: c «= a 1 r + a 2 

With the values of c at r = 1.25 and r = 2.1875 the constants 
a 1 and a 2 are found as follows: 
0.281 = 1 .25 a 1 + a 2 

0.217 = 2.1875 a 1 + a 2 

0.064 = -0.9375 a. 
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-0.068 

-0.068 x 1 .25 + a 2 
0.366 

-0.068 r + 0.366 

Suppose we have the foot of the blade at r =0.2 R (at smal- 
ler radii the contribution of the extracted energy is neg- 
lectable, the more so in relation to the material requir- 
ed for the blades), then we can calculate the chords at the 
foot and at the tip of the blade. 

c foot = -°* 068 x °-5 + 0.366 = 0.332 
c + - = -0.068 x 2.5 + 0.366 = 0.196 



In order to see which deviation is introduced by lineari- 
zation of the blade at the various cross sections, the com- 
puted values for the chords are given in the following table. 



cross 


r[m] 


c [mm] 


c [mm] 


c[mm] 


deviation 


section 




exact 


first 


second 


[mm] 


number 






linearization 


linearization 




1 


0.5000 


252 


332 


250 


2 


2 


0.6250 


273 


324 


(264) 


9 


3 


0.9550 


293 


301 


300 


7 


4 


1 .2500 


281 


281 


(281) 





5 


1.5625 


260 


260 


(261) 


1 


6 


1 .9100 


235 


236 


(238) 


3 


7 


2.1875 


217 


217 


(222) 


5 


8 


2.5000 


198 


196 


200 


2 



Prom thi3 table it can be seen that considerable deviations 
occur at the two bottom sections of the blade (cross section 
1 and 2). 

Therefore a second linearization will be useful for those two 
sections, the more so in relation to the material requirements 
for- the blades. At the same time rounded figures can be chosen 
for the relevant measurements of the blade. 



— a 1 = 
0.281 = 

— *2= 
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So the first linearization provides us with c = 200 mm for 
the tip and c = 300 mm for cross section number 3, while the 
second linearization provides us with c = 250 mm for the foot 
of the blade. 

The actual deviations then existing are shown in the last 
column of the above table. 

Linearization of the blade angles /3 (or the blade twist) can 
be done in a similar way. However, here we first have to de- 
cide how the blades will be fixed to the rotor frame. In ca- 
se the blades are fixed to the rotor frame a 1 ; two positions 
then only one linearization is possible since a straight line 
is defined by 2 points. 

Suppose the blades are supported in cross section number 3 

and 6 (which is about at 1/4 and 3/4 of the blade length) then 

the angles /3 can be found as follows: 

ft can be written in the following way /3 = a^r + a ^ 

With the values of (b at r = 0.955 (cross section 3) and at 

r = 1.910 (cross section 6) the constants a^ and are found 

as follows: 



32.2 




a 3 x 0.955 


+ a 4 


19.1 




a^ x 1 .910 


+ a 4 


13.1 




-0.955 a 3 




a 3 




-13.72 




32.2 




-13.72 x 0. 


955 + a 4 


a 4 




45.3 




So (b 

t 




-13.72 r + 


45.3 



The computed values for /3 at the various cross sections are 
presented in the following table. 
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cross 

section 

number 


rim J 


exact 


Q r Oi 

linearized 


deviation 

r Oi 
[ ] 


1 


0.5000 


42.6 


38.4 


4.2 


2 


0.6250 


39.3 


36.7 


2.6 


*< 




^2.2 


^2.2 


o 


4 


1.2500 


27.0 


28.2 


1.2 


5 


1.5625 


22.8 


23.9 


1.1 


6 


i .9100 


19.1 


19.1 





7 


2.1875 


16.8 


15.3 


1.5 


8 


2.5000 


14.7 


11 .0 


3.7 



The result of the linearization is visualized in the 
following graphs where blade form and blade twist are 
shown. 



[ 



rotor centre 



+ 



I 

cross section: 1 



i 

4 



I 

T 
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5.3 Calculation of the radius of the curve of the blade 

The radius of the curve of the blade is depending on the chord 
of the blade. For a 10$ arched steel sheet (f/c =0.10) the 
relation of radius and chord can be derived as follows. 




From this figure we will see 
that in the right-angled tri- 
angle the following relation 
is valid 

R 2 = (R - 0.1 c) 2 + (1/2c) 2 



Worked out it will give us 



R = 


1.3 c 




(5.8) 


For 


a 7$ arched steel 


sheet 


the 


relation will be 




R = 


1 .8 c 




(5.9) 



In our calculation example we assumed the blade to be support- 
ed in cross section 3 and 6. The corresponding chords are 300 
mm resp. 238 mm. So the radii of the curves of the blade at 
those positions and of course of the blade supports as well 
have to be: 

in cross section 3 : R = 1.3 x 300 - 390 mm and 
in cross section 6 : R = 1.3 x 238 = 310 mm. 
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5.4 Measurements of the blade . 

In our calculation example we found the chord of the blade. 
This is, however, not the cutting measurement of the blade. 
Therefore we have to calculate the length of the segmental 
arc with the computed chords. 

For a 10^ curved blade we can calculate this as follows: 




This correction of 2.6$ is to be neglected. 
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6. DESIGN AND STRENGTH ASPECTS OP THE WINDMILL 
6.1 Strength calculations . 

6.1.1 Strength calculations of the rotor of the 
12 FU 500 windmill. 

On the blades are acting the centrifugal forces and the 
thrust forces of the wind. For the calculation of the 
centrifugal forces we assume a wind velocity of 20 m/s 
and an unloaded windmill. This means that the tip-speed 
ratio will be about 1,65 times more than in case the 
windmill is matched to a proper load. 

The angular speed will be: 

AB k s 1.65 x 2 x 20 . 26 rad/sec 

The formula for the centrifugal force is 

F c = S~ or F c = m si 2 r 

The mass of one blade is approx. 4 kg and the centre of 
gravity of the blade is situated at a radius of approx. 
1 .5 m. 

P c - m.r?.r = 4 x 26 2 x 1.5 ■ 4056 N 

To calculate the number of bolts one has to fit the blade 
with, we first have to calculate the maximum load which 
can be taken by one 1/4" bolt. 

The area of a 1/4" bolt is -5- 4.8 2 = 18 mm 2 

The allowable unit stress for steel is 120 N/mm 2 

So one bolt can take 18 x 120 = 2160 N 

The minimum required number of bolts has to be jffifl = 1.6. 

In fact the complete centrifugal load of the blade is taken 
up by the inner ring, since this one is much more rigid than 
the outer ring. 

Although 2 bolts should have been sufficient, from the 
strength point of view, there are 3 bolts applied in order 
to maintain a proper curve of the blade. For the outer ring 
2 bolts have been applied in order to support the blade 
properly..., , , .. .. . ............... ' 
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6.1.2 Strength calculations of the tower . 

The major forces acting on the tower are the thrust force 
of the rotor, the weight of the total head construction, 
tail and rotor and the pumprod forces. 

To calculate the thrust force created by the wind we apply 
the following formula: 

P p = | x 1/2 p V 2 A (6.1) 
Fp = thrust force (Newton) 



8 
9 



= Thrust coefficient 

^ = Density of air (1.2 kg/m^) 

V = Wind velocity (m/s) 

A = Swept area (m ) 

Although the windmill is put automatically in secured 

position if the windspeed exceeds 10 m/sec, we will make 

the calculation for a windspeed of 20 m/sec. So even 

in case the security mechanism fails, the strength 

of the system (and the structure) will be sufficient to 

withstand the forces acting on it in a windspeed of 20 m/s. 

Note that for normal operation, in a windspeed of 10 m/s, 
a safety factor of 

20 2 

— — = 4 is obtained. 
10 2 



For a rotor with a diameter of 5 m the thrust force will be: 
P p = | x 1/2 x 1,2 x 20 2 x J|- 5 2 = 4189 N 



The weight of the total head construction, tail and rotor of 
the 12 PU 500 windmill is approx. 300 kg. This results in a 
vertical force of 3000 N , acting on the tower. 
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Pumprod forces are assumed not to exceed 3000 N during the 
upward movement of the piston. 

So the total vertical force acting on the tower amounts to 
6000 N. 

In the following diagram these forces are situated at the 
rotor centre. 

The thrust force causes also a bending moment acting on the 
tower which is also shown in the same diagram. 

uic't^KV of h*audLcon%VruLc\c'«oT\ t roVor «ndL 
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With help of Cremona diagrams the member forces of the 
lattice tower have been derived on the fore going pages. 
These member forces are derived for each seperate load 
acting on the system. The final resulting member forces are 
as follows: (in Newtons) 



No 


Tension 


Thrust 


1 


6550 




2 


6550 




3 


17825 




4 




23975 


5 




27050 


6 




27050 


7 


9125 




8 


6525 




9 






10 






11 


3175 




12 


10775 




13 




2312 



It can be seen that the maximum thrust force occurs in the 
members 5 and 6 (and in 1 and 2 if the wind direction is 
opposite). Because it is a four-leg tower this force is 
distributed among 2 legs. 
So the thrust load on one leg is 13525 N. 



The unit stress in compression of the angle irons of the 
tower legs can be calculated as follows: 

For angle iron 40 x 40 x 5 mm, we can find from handbooks 
the following data: 

cross section area A = 3*79 cm 2 
Moment of inertia I v «= 2.26 cm^ 
Radius of gyration i y = 0.77 cm. 

The maximum buckling length is 2 m (l = 200 cm) and the 
members can be considered as fixed in between rigid joints. 
So the slenderness ratio will become 

X - ILL - TOO - no 
A i7 "0777 130 

For normal construction steel the allowable unit stress in 

compression for this slenderness ratio of 130 will be 

(according handbooks) C7^ = 3500 N/cm 2 . 
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The true stress is given by: 



buckling force 
cross section area 



= 3569 N/cm 2 . 



This is approx. the maximum allowable. 

On the diagonals of the tower (members 7 to 12) only tension 
forces occur with a maximum of 10775/2 « 5388 N for member 12 
(and 11 for winds from opposite direction). 
These members consists of angle iron 30 x 30 x 4 with a 
cross section area A of 2.27 cm 2 . 

The true stress set up in these angle irons is given by: 



The allowable tension stress for construction steel is 



It can also be noticed that the members 9 and 10 are not 
carrying any load. They can, however, not be omited because 
they halve the buckling length of the members 1,2 and 5 t 6. 

We may conclude that this tower design can withstand even 
those forces which would occur with a windmill in working 
position at a wind speed of 20 m/s, whereas the normal limit 
of operation is at 10 to 12 m/s. 



C= 



tension force 



cross section area 



4S-= 2374 N/cm 2 . 
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6.1.3 Bearings . 

Selection of the proper bearing for a given application 
may be a difficult problem, but some general principles 
can be set down. 

Ball bearings are usually the most economical choice for 
the smaller sizes and lighter loads; roller bearings are 
usually more satisfactory for the larger sizes and heavier 
loads. Roller bearings have considerably less sensitivity 
to shock loading than ball bearings. Most bearings lose 
capacity if misaligned any measurable amount, but self- 
aligning bearings can take 2 or 3 degrees of misalignment 
without loss in capacity. 

A major influence on bearing choice is the direction and 
magnitude of the load. Self-aligning ball bearings are 
designed primarily to carry light and moderate radial 
loads, but can carry some thrust load in either direction. 
Deep-groove ball bearings are designed to carry light, 
medium or heavy radial or combined loads. These bearings 
can carry substantial thrust loads. 

In the following table the characteristics of a single 

groove bearing and a self-aligning bearing with a 50 mm 

bore (medium series, easily available) are given (SFK-bearings ). 



50 mm bore. 


Single groove 


Self-aligning 


Dyn. ref capacity C 
Thrust factor Y 
Rotation factor X 


47000 

1-1.5 
1 


39000 
3.2 
1 



The thrust force acting on the rotor of the 12 PU 500 
windmill at a windspeed of 10 m/s is calculated with 
(6.1) as Pas 1047 N. 

We assume that this thrust force is completely taken by 
the front bearing. 
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rear 
bearing 



The radial load acting on this bearing is caused by the 
weight of the rotor and is calculated as follows: 

rotorweight 
P A 161 

front 16, 
bearing ^ 

PB 

The rotorweight amounts to 161 kg and the centre of gravity 

is located at approx. 0.1 m from the front bearing. 

The distance between the two bearings is 0.5 m. 

According the conditions of equilibrium, the reaction forces 

can be calculated from: 

£ M = o 

With regard to A: 1610 x 0.6 - P fi x 0.1 = 

P B = 1930 N 

and from ^H = o F A + 1610 - 1930 = 

P A e 320 N 

Now the equivalent thrust load, defined as 

P = X.P r + Y.P a (6.2) 

can be calculated for the front bearing as 

P = 1 x 1930 + 1.25 x 1047 = 3239 N. 

At a windspeed of 10 m/s the revolution speed of the 

rotor will be approx. 110 RPM. , ^ 

So the RPM factor becomes : f Q = \J - \J ^1 - o.67 

and f - f "« C - 0.67 x 47000 — q ho 
8,1(1 *h ~ p" " = 9.72 

This results in a lifetime Lh according to: 

% - 500 fl (6.4) 

Ljj = 500 x 9.72 3 s 459165 hours. 

With an average of 5000 running hours per year the life of 
the bearing will be about 92 years. This is, however, only 
valid if the bearing is running under proper conditions, 
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well greased in dust free housings, etc. Under field con- 
ditions the life will be reduced depending upon the degree 
of maintenance, proper alignment, etc. 

For the calculation of the rear bearing, we have two 
radial loads acting on the bearing. One is caused by the 
weight of the rotor and has been calculated earlier as 320 N. 
The other is caused by the pumprod forces. 

The magnitude of these pumprod forces is not exactly known, 

but is assumed to be 5000 N during the upward movement of 

the piston. During the downward movement of the piston the 

direction of the pumprod forces is opposite, but also limited. 

So the maximum load acting on the rear bearing is: 

5000 - 320 = 4680 N. 

The equivalent thrustload is now 

p = x .F r + YFa = 4680 N 



So Ljj = 500 f h J = 500 x 6.73 J = 152410 hours, 
which means a life of 30 years. 

For the self-maligning bearing these values are for the 
rotor side bearing calculated as 1 2 years and for the 
rear bearing as 17 years. 

It may be concluded from these calculation examples, 
that both types of bearings may serve the purpose. 
However, in case single groove ball bearings are used, 
proper alignment is most essential for the life of the 
bearings (and a straight shaft as well of course). 



and f , 



f n .C _ 0.67 x 47000 



= 6.73 
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6.2 Flywheel energy of the rotor . 

A flywheel is primarily a rotating energy reservoir 
performing two distinct classes of service. First, the 
flywheel absorbs energy from a power source (luring the 
greater portion of the operating cycle and delivers a 
large amount of energy as useful work in a very short 
portion of the cycle. This class includes punch presses 
and shears. In the second class, including steam engines, 
internal-combustion engines, compressors and reciprocating 
pumps, the flywheel smooths out the speed fluctuations 
caused by non-uniform flow of power from the piston 
during each energy cycle. 

A windmill rotor is also to be included in the second 
class where the rotor is also acting as a flywheel in 
combination with the reciprocating piston pump. 

The kinetic energy E (J) in a rotating flywheel is 

E = f W U 2 

Where in W = Weight of wheel (kg) 

U = Velocity at radius of gyration (m/s) 

If the velocity changes from U-j to U 2 the energy fluctuation 
is: 

A E = i W (U* - u|) 

We will make a calculation example for the 12 PU 500 windmill 
at a wind velocity of v = 6 m/s . 

The weight of the complete rotor is 161 kg. The radius of 
gyration is assumed as 1.25 meter. The tangential velocity 
at r = 1.25 m in a wind velocity of v s 6 m/s and with a 
tipspeed ratio of X*: 2 is: 

D = | . X. v = x 2 x 6 = 6 m/s. 

So E = |WU 2 = |x 161 x 6 2 = 2900 J. 

The energy fluctuation can be determined by calculating 
the generated energy as follows: 
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The torque generated by the windmill at a wind speed v 
is: (relation 4.10 and 4.11) 

T = ipv 2 irR 3 . C t and C t = ^ 

For \ = 2 — ► C t = P-£2§ = o.19 

So T = i x 1 .2 x 6 2 x tt x 2.5 3 x 0.19 = 201 Nm. 

For one revolution this presents an amount of energy of 
201 x 2TT = 1263 J 

Summarising this we conclude that the energy stored in the 
flywheel (= rotor) is 2900 J and that the energy generated 
during one revolution is 1263 J t which is the maximum the 
pump can "ask", (at a wind speed of 6 m/s) at an average. 

In the following graph this is made visable. 
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6.3 ^owerheight . 

The cost price per kwh (Cg) generated by a windmill is 
given by the following relation 

n _ Total initial investment + rec.cost + mainten ance cost 
Total energy generated during the same period. 

For the various existing windmills it has been discovered 
that a minimum for Cg is obtained for 1.2^jj < 1.7 

H = Height of tower 

D = Rotor diameter 



The actual costprice per kWh is naturally depending on the 
actual wind pattern in a particular area. 



If we make the calculation for a 12 PU 500 windmill in 
Ghazipur area, we will get the following result, with 
the following assumptions made. 

Cost price of windmill Rs. 7000 

Lifetime of windmill 15 years 

Yearly maintenance cost Rs. 200/year 

Yearly pumped quantity of water 20000 m3 

Over a total head of 6 metres 
Rate of interest 10# (i=0.10) 



We have to calculate with an annuity factor (a) of 

a = o.iox (1 +Q.1Q) 15 = . 131474 



= i * + i) n = a = 0.10 x (1 + 0.10) 1 5 
(1 + i) n - 1 ~ (1 + 0.10) 15 - 1 



The total energy generated can be calculated as follows: 
(in kwh) 

E = ?-g- H -Q g 1000 x 10 x 6 x 15 x 20000 _ 500Q km 

1000 x 3600 1000 x 3600 

If we apply the formula for the cost price calculation per 
kwh it will give us: 

n _ 0.13147 x 7000 x 15 + 15 x 200 

5000 

1 ^ 804 * ? 000 . IMP! , Rs. 3.36 / kwh 
' — - *» nf>n - - * nr>r> 
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It must be stressed that this kWh-price (as calculated for 
the windmill ) should not be compared directly with kWh- 
prices of electricity. All losses of the system (pump, 
rotor losses) have been taken into account in calculating 
this price. With 1 kWh electrical energy an electric motor- 
pump unit must be driven which has an efficiency of less 
than 5096. The real kWh-price is thus at least twice as high. 
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7. WIND ANALYSES 

To investigate what type of windmill will be the most 
suitable one in particular area and to estimate the eco- 
nomic benefits which can be derived by application of 
windmills for waterpumping, analyses of the wind 
potential in the area is the main instrument. This chap- 
ter will deal with the way, how to analyse wind data for 
this purpose. 

7.1 Availability of data 

Often wind data are registrated at meteo-stations and at 
airports. At most meteo-stations, however, only daily or 
monthly averages are collected and published. These long- 
time averages are only of limited value for our purpose 
as we will see later. At airports mostly hourly wind- 
speed data are registrated during day time, when the 
airport is in operation. Only at big international air- 
ports continuous registration is practised. Sometimes it 
is hard to obtain these data, since often these data are 
not published. As we will see further on in this chapter 
it is essential, to have data from a number of years 
to get an impression about the variability of the wind 
in the area. 

7.2 Reliability of the data 

Wind is one of the most fluctuating meteorological fea- 
tures and influenced by all types of disturbances topo- 
graphical hinderances as well as air disturbances. 

Since it takes often too much time to collect wind data 
in the place where the windmills will be installed, 
one often has to rely on data from a meteo-station or 
airport quite far away. Roughly it can be said that when 
no major obstacles occur -(mountains) and both places be- 
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long to the same climatological area the geographical 
variation in the wind speed will not be very large. 
At least the data obtained from the nearest station 
will give an indication about the wind potential in the 
area. A better impression can be obtained with short 
period on-site measurements with which the degree of 
correlation between the site and a nearby established 
wind speed recording stations wind history can be deter- 
mined. 

The extent to which site wind patterns reflect regional 

patterns is the degree of correlation. 

Traditional site evaluations show that in case a high 

degree of correlation exists, then on-site wind speeds 

may be projected into the future based on the nearby 

station's past windspeed records corrected by a constant 

factor. 

The reliability of the data, depends largely on the way 
of observation. Meteo-stations are often situated near 
or close by cities and surrounded by buildings or trees. 
Moreover the wind speed is often observed at 2m height. 
The figures obtained in this way, will be far below the 
wind speed in an open area, where windmills will operate. 
Airports are always situated in open areas and the wind- 
measuring equipment is often placed on top of the air- 
traffic control building. Except for the difference in 
height of observation for which a correction has to be 
made, these data are proper for further analysis. 

It is advisable not only to collect the data, when avail- 
able, but also to see how the data are collected. Continu- 
ous observations by a recorder are much more reliable than 
momentaneous observations say once in every hour of the 
wind speed. The latter is no indication of the average 
hourly wind speed at all. However, if every hour a wind run 
indicator reading is taken, this provides also an hourly 
average . 
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7.3 Annual Wind Pattern 

The annu .1 wind pattern can provide some information 
about the variability of the average wind speed during 
a year. In this data the first indication can be found, 
whether the wind pattern coincides with the annual pat- 
tern of crop irrigation requirements. 

Graphical presentation of the average wind pattern of 
Babatpur airport (Varanasi). 




1 • 
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Here it can be seen that the period with highest wind 
speed coincides with the summer season, when a lot of 
irrigation water is required to grow a crop. 

7.4 Diurnal Wind Pattern 

The diurnal wind pattern is the wind speed distribution 
during the day. On the next page the diurnal pattern for 
every month of the year can be found, derived from data 
from Babatpur Airport (Varanasi). From this type of 
curves it can be found on the average at what periods of 
the day, the windmill will run and when the peak output 
can be expected. This is an indication whether the wind- 
mill output coincides with the daily requirement and if 
not what capacity the storage tank has to have to store 
the water pumped during night time for example. 

Since the annual windspeed pattern as well as the diurnal 
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pattern are average figures it is not possible to derive 
from them the capacity of a windmill pump system. 
To estimate the real amount of energy which can be 
extracted from the wind, the frequency of occurrence of 
periods with a certain wind speed has to be counted. 
With this we will deal in detail below. However, we will 
first study the maximum wind speed figures. 

7.5 Maximum wind speed 

The maximum wind speeds are important for the strength 
calculation of the windmill structure and for the design 
of the security mechanism. If cyclones or other heavy 
storms occur in the area, proper protective measures 
have to be taken to minimize the chance of damage to the 
windmills . 

Designing a structure which can withstand all possible 
wind speeds will be impossible and also costly. Again it 
is stated here that hourly average maximum wind speed 
figures give no indication of the real peak gusts 
which may occur during a few seconds, but can be very 
destructive. 

7.6 Windspeed frequency distribution 

The wind speed frequency distribution is the most im- 
portant wind characteristic for optimal windmill design 
and for windmill output predictions and consequently for 
windmill economics. 

Hourly average wind speeds are the minimum requirement 
for this approach. On the next page one example has been 
given in what way the data have to be processed to derive 
this frequency distribution. 
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RECORDED WIND SPEED 



Yoor : 1978 
Month: May 
Day : 17 th 




First the hourly average is calculated out of the recorder 
readings, from the observation station. 

These hourly averages are collected in a table as shown 
below. From this table which contains all the hourly 
averages of one month, also the hourly mean, daily mean 
and monthly mean can be derived. 



12 
13 
14 
IS 
16 
17 

31 
hourly 



r 1978 
1 2 


3 
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5 
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7 


8 
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9 10 


11 
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16 


17 
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19 


20 


21 


22 


23 


doily 
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1 


1 


1 

1 




1 
1 


1 
1 


1 


1 


1 


1 


1 


i 


1 


1 


1 


1 


1 


1 


1 






1 


1 


1 


! 1 


3.2 


3.9 


3.9 


3.0 


2.7 


2.6 


2.9 


2.9 


2.2 


1.1 


1.3 


1.0 


1.8 


2.2 


1.3 


0.9 


1.7 


1.7 


1.9 


2.5 


2.9 


2.4 


3.0 


3.9 


2.4 


4.2 


4.3 


5.0 


3.9 


3.6 


4.2 


3.6 


3.2 


2.9 


2.1 


2.4 


3.2 


3.2 


2.9 


2.7 


2.0 


2.7 


3.6 


2.0 


2.0 


3.0 


3.7 


2.3 


1.9 


3.1 


2.1 


1.4 


0.6 


0.7 


1.2 


1.4 




1.4 


1.4 


2.2 


2.3 


4.0 


4.3 


4.5 


4.5 


4.2 


4.0 


3.6 


3.0 


2.4 


2.2 


2.3 


2.4 


2.9 


2.9 


2.4 


2.2 


1.1 


0.9 


0.4 


1.0 


1.1 


2.0 


2.9 


2.4 


1.4 


1.4 


1.4 


2.2 


1.9 


1.9 


2.1 


1.4 


1.0 


0.7 


1.2 


2.9 


3.0 


2.9 


,, 


3.0 


2.5 


2.2 


2.6 


2.9 


2.7 


4.3 


4.8 


5.0 


5.0 


5.0 


5.3 


5.2 


5.2 


5.0 


5.2 


5.1 


5.2 


4.3 


2.8 


3.5 


4.0 


4.3 


4.3 


4.1 


3.9 


3.2 


4.3 


4.3 


3.2 


3.2 


5.0 


5.3 


6.0 


6.4 


6.4 


6.3 


6.2 


6.4 


6.3 


6.2 


6.4 


6.0 


5.0 


4.3 


3.7 


3.5 


2.7 


2.2 


4.9 



































































































[3.1 


2.9J2.9 


2.6|2.3 


2.5 


3.4 


3.5 


3.4 


3.2 


3.1 


3.5 


3.7 


3.9 


3.6 


3.4J3.7 


3.6 


2.9 


2.5 


2.8 


3.1 


3.0 


3.0 


3.1 


monthly 
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Secondly the occurrence of wind speeds in each range of 
wind velocities with the intervals of 1 m/s are counted 
and collected as shown on the following table. 



Month : May 1978 



wfndipMd 
tang* m/» 


fraquMcta 


In 

haun 


in 

{MrcantagM 


12 - co 








11 - 11.9 








10- 10.9 








9 - 9.9 








8 - 6.9 








7 - 7.9 








6- 6.9 


tt it- u& 
fill 4WT 


10 


7 


5- 5.9 


m m mi 


14 


10 


4- 4.9 


w mm w 


18 


13 


3 - 3.9 


m m mwt m i 


26 


18 


2 - 2.9 




44 


31 


1 - 1.9 
0- 0.9 


40 m w **■ »t 


25 


17 




6 


4 




total 


143 


100 



!Ehis frequency distribution is visualized in the follow- 
ing Graph, a so-called histogram. 



Windspeed frequency distribution 



Li 



p.: 

r 



m 



Hi 



-i — i 

> m/s 



-r- 
8 



6 » 8 3 

windspeed - 



10 



5 3 



When this has been completed for every month of the year, 
these data need special observation to find the most op- 
timal windmill pump system. 

The design wind speed of the windmill is determined by the 
frequency distribution of the wind speed and the output 
characteristic of the windmill pump system. 

7.7 Determining the design wind speed 

The performance of the 12 PU windmills with a reciprocating 
single acting piston pump can be approximated by a graph 
as indicated below. /p ^y*^ 

Ideal output 



Real output 12 PU 
Approach 



0o 




■O v 



For V = Yd may write: 



or: 



P = P 

windmill water 

7tot J? V d 3 A = P. « " f 

,/ytot jf v a ij_ 
? » « » 



(7.1) 
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For V>V d we may write: 

(a straight line through the origin and the intersection 
of V d - ft) 

. ^tot^fVd 2 A 
Pw* H 



= 



. V 



(7.2) 



For the 12 PD windmill the overall efficiency (7 tot ) at 
the design windspeed is 0.3 
With f> w = 1000 kg/m3 

g = 9.8 m/s2 

f = 1.2 kg/m3 



we may write for (7.2) : 

= 1.84 » 10~ 5 ^ x V d 2 x V [m 3 /s] 

The value of V d will be optimal for the maximum of 



(7.3) 



Q = f"0dt 
j ti 

Substitution of 7.3 in 7.4 gives: 
Q 



(7.4) 



" 5 H K V d 2 " v dt 



= J^ 2 1.84*10' 

A, 2 

Q = 1.84«10~ 5 « £ « V d 2 J v 
In other words Q is maximal for V 



dt 

2 



(7.5) 

dt = maximal 
(7.6) 



With the wind speed frequency distribution it is easy to 
find this optimum, if we write (7.6) in the following way: 



'1 



2s* 

dt ^ v d < Vi ti 
*1 



(7.7) 



Vd = design wind speed 

Vi = interval middle wind velocity 
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With the following example this procedure will be 
explained : (same example as before) 





V 


Vd 


Vi 


ti 


Vi ti 


t 2 
> Vi ti 

^1 




2 .J 2 
V|^ Vi ti 

*1 


o 


- 1 





0.5 


6 


3 


464.5 








1 


- 2 


1 


1 .5 


25 


37.5 


461 .5 


1 


461 .5 


2 


- "* 


2 


2.5 


44 


110 


424 


4 


1696 


j 


- 4 


-> 




26 


91 


314 


9 


2826 


A 


- 5 


4 

"T 


4.5 


18 


81 


223 


16 


3568-4 


5 


- 6 


5 


5.5 


14 


77 


142 


25 


3550 


6 


- 7 


6 


6.5 


10 


65 


65 


36 


2340 


7 


- 8 


7 


7.5 












8 


- 9 


8 


8.5 












9 


-10 


9 


9.5 












10 


-11 


10 


10.5 













Prom the last column can be seen that the maximum value 
for Q is reached for Vd of 4 m/s; or better that the design 
wind speed should preferably lie between 4 and 5 m/s in 
order to achieve a maximum output for the given frequency 
distribution. 

It is stressed here that the optimal wind speed over the 
year can vary considerably. The design speed should be 
chosen in such a way that the maximum obtainable output 
will be reached during the critical months. 
To obtain also optimal outputs during months with lower 
wind speeds the stroke of the pump can be reduced by 
reducing the acting crank arm. 

Experiences reveal that it is more important for the 
introduction of the windmill that the windmill works 
and pumps at least something, than when it only works 
during a few hours of the week. 
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If a regular water supply is regarded to be more important 
than maximizing the water yields, the design wind speed 
should be chosen at a lower value, for example from 3 to 
4 m/s. 

In case of the example used here with a chosen design 
speed of 4.5 m/s the windmill will run only 23# of the time. 
With a design speed of 3.5 m/s the windmill will run 38$ 
of the time with an approx. 20$ reduced output. 

The frequency of wind speeds of 10 m/s and more should 
not be taken into account in determining the optimal 
design speed since the windmill will be in secured 
position at those wind speeds. 

Concluding we may say that the maximum obtainable output 
is not the only main objective. 

Also the running percentage has to be considered in 
connection with the required regularity of water supply 
and for psychological acceptance of the windmill. 
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8. WINDMILL OUTPUT CALCULATIONS 



If the input-oot put relation of a certain windmill is 
known, the total quantity of water which can be pumped 
during a certain period can be calculated with the 
help of the wind speed frequency distributions. This 
is, however, only valid for one combination of pump- 
size, total head and design wind speed. 

The theoretical output can be calculated with (7.5) 
and (7.7) for a windmill of the 12 PU serie: 

t 2 

Q = 1.84 . 10~ 5 $ V d 2 ^ Vi.ti 

For ti in hours instead of seconds this relation becomes: 

*2 



Q = 6.624 . 10"* • v d ^ Vi ti (8.1) 

t 2 
t 



The tf>rm V d 2 ^ Vi ti is calculated in the table of chapter 7.7. 



1 

as 3568 for a Vd of 4 m/s. 

So the theoretical output during the period observed with 
a 12 PU 500 windmill and a total head of 6 m will be: 

2 

Q = 6.624 * 10~ 2 * TT * 2 ' 5 * ^368 

5 

Q = 928 m 3 

In practice however, this value will be always more or less lower 
due to restrictions in the number of different diameters 
of the standardized pumps. One can make a choice out of 
3% 4 M , 5 W , 6" and 8 M diameter cylinders. A windmill with 
an oversized pump results in a higher design speed and 
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consequently in a lower output (and fewer running hours). 
On the other hand a windmill with an undersized pump 
will result in a lower design speed and consequently in 
a lower output (but more running hours). 
In chapter 9 determining the diameter, of the pump will 
be dealt with. 

Also variation in ground water table (or other water 
source) will irfluence the design wind speed and conse- 
quently the output. 

Often there is a positive correlation between a re- 
ceding water table in summer season and an increasing 
wind speed* The effect of increasing wind speed is quite 
often even sufficient to compensate not only the in- 
creased lifting head, but also to meet the higher water 
demand during summer season. 

These three effects have to be taken into account in 
determining the most critical month (s). 
Comparison of theoretically calculated outputs and the 
actual windmill outputs revealed a difference of about 
30$. The actual output being 70$ of the calculated output. 
This ratio is called the working efficiency of the wind- 
mill and is estimated to be about 70$. 

It is due to the fact that the windmill will work a less 
number of hours than theoretically possible due to mainte- 
nance, breakdowns or being left in secured position (semi- 
automatic security mechanism) (v. Vilsteren 1979) and 
due to fluctuations in wind speed and changing wind 
directions. 

Finally it will be clear that wind data of one year can 
not provide sufficient reliable information concerning 
future wind energy potential. At least 5 to 10 years of 
data have to be worked out to find reliable data by means 
of averages and probability analyses. 

Although collecting 5 to 10 years of on-site wind data, 
while extremely desirable in terms of providing high 

quality data, would most likely be prohibitively expen- 

- J - ' -~ - - : - - - - * ■ ---- ■ ■ ■■ 
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8.1 Correction factor for difference in height 

The height of the rotor shaft above ground level is seldom 
the same as the height at which the windspeeds are recorded. 
Hence a correction factor has to be applied for the wind speeds 
"seen" by the windmill rotor since the actual wind speed is 
depending on the height at which it is considered. 

The friction between the lowest layer of the atmosphere and 
the earth is largely determined by the roughness of the earth's 
surface. A smooth surface like that of a beach will have little 
effect upon the wind speeds in the layer above, but a rough 
surface, like a city, has a strong effect. 

A secondary effect is the stability of the atmosphere. 
One can imagine that the turbulent exchange is much higher 
in an unstable atmosphere, with the result that higher wind 
speeds in upper layers more easily can reach the surface. 
This effect becomes smaller, however, when the wind speeds 
are higher. Because we are mostly interested in higher wind 
speeds, for wind energy purposes we will neglect this effect. 

Generally, meteorologists have adopted a logarithmic profile 
to describe the increase in average wind speed with increasing 
height above the roughness elements. 

Usually one wishes to determine the wind speed at a height H 2 
when the wind speed at a height is given. 
The expression therefore is: 

H2 

V H2 z o 



v m mir (8 - 1) 

z 

z is the roughness height. 

This assumes, however, that the roughness heights at both places 
are the same. 

When the roughness heights differ, then we need a reference 
height where we can assume that the wind speed is not much 
affected by the roughness. This reference height is often 
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determined as 60 m (above the roughness elements). 



The formula then becomes: 



v m §2- x m P- 

V H2 = Z Q1 Z Q2 

V H1 in SSL x in S- 

z o 2 z oi 



(8.2) 



When the roughness heights are the same (8.2) reduces to 
(8.1). 

The roughness heights are difficult to determine. 

A number of values of z for typical terrain types is given 

in the following table. 



Terrain description of area within several 
kilometres upwind of site 



z m 



smooth surface, open sea, desert (flat) 
low grass, flat open terrain 
high grass, airports (runway area) 
farmland, isolated trees, small crops 
farmland, many hedges, tall crops 
many trees, hedges, few buildings 
small towns, suburbs 



0.001 

0.008 

0.02 

0.03 

0.08 

0.2 

0.6 



The international accepted standard height for taking wind 
data is 10 m. 



We will make a calculation example for the following situation, 
wind speed recordings have been obtained from an airport 
with a cup anemometer at 17.2 m above ground level (H2 - 17.2 m 
In case of applying a windmill with a rotor height of H1 = 7 m 
above groundlevel situated in farmland vdth tall crops, we 
may take the following roughness heights according the table 
for the airport z O2 =0.02 m and for the windmill site z 01 =0.08 m 

The correction factor will become: 



m 60. x ln |JL 

V H1 z 02 Z 01 



Z 01 z 02 



1r , 60 v , 7 

= ln 60 ln 17.2 = °' 8 ° 



So the wind speed at rotor height is 80 $ of that at measuring 
height . 



It can also be seen here that increasing the height of the 
windmill (tower) will have a positive effect on the power 
which will be generated, but a negative one with regard to 
the cost price of the tower as we have seen already in 
chapter 6.3. 
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9. DETERMINING POMP DIAMETER 



In order to determine the diameter of the pump (cylinder) 
one has to know the total head over which the pump has to 
lift the water and the design wind speed, which resulted 
from the wind speed analyses (see chapter 7.7). 
For V= Vd we may write (see 7.7) 

P = P 

windmill water 

ifV d 3 TTR 2 . 7tot = p w .g.H.V fa . n . % (9.1) 

where: p = air density = 1.2 kg/m^ 

V d = design speed [m/s] 

R = Rotor radius [m] 

p w = density of water = 1000 kg/m 3 

g = acceleration of gravity =9.8 m/s 2 

H = lifting head [m] 

n = number of revolutions [rev/s] 

V g = stroke volume jm 3 ] 

"?tot ~- tota l efficiency 

"7 V ~ volumetric efficiency 

For the number of revolutions of the pump we may write 
at V = V d : 

» = (9.2) 
2TTR 

and for stroke volume of the pump 

V s = ^TTd 2 x s (9.3) 

where: s = stroke of the pump [m] 

d = diameter of the cylinder [m] 
Substitution of (9.2) and (9.3) in (9.1) results in: 

d 2 - a £ JL R 3 • toot V 2 , Q „ , 

d " 4 F w - T ' IT 7v~ ' A (9 ' 4) 
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For the 12 PU series we write for: 

X = 2 
"7tot = °'30 



So that (9.4) becomes: 



Vh.s 



d = 0.0165Vd\/ii — (9.5) 

H.s 



and for the 12 PU 500 with R = 2.5 m 
d _ 0.065 Vd 



(9.6) 



For the example which was taken in 7.7 we found a design 
wind speed of 4 - 5 m/s, say 4.5 m/s. 

So for an elevation head of 6 m (which was taken in the 
example of 8) we find for the diameter of the pump, for 
a stroke of 0.24 m: 

d . 0. 065 x 4.5 . 0>24 m 
\J6 x 0.24 

If, however, the design speed is chosen lower, say 3 m/s, 
for reasons mentioned in 7.7 the diameter will become: 

d = °- 65 *3 = 0.163m 
\j6 x 0.24' 

This results obviously in a 6" ••standard" pump, with a 
herewith belonging design wind speed of 

0.15 = °' Q 65 x Vd => V d = 2,8 m/s 

x 024 
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If this pump (6 M ) is coupled to the windmill then the 
expected output will follow from the table of 7.7. 

Bie term V d 2 ^ 2 Vi ti is for 3 m/s: 2826, so with 

expression 8.1 the output during the observed period in 
the example will be: 612 m^. 

Biis is 34$ less than the theoretical optimum, but the 
number of running hours is increased from 239S to 48$ of 
the total time and so a more regular delivery is obtained. 
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10. MEASURING OF WINDMILL PUMP PERFORMANCE 

After the windmill has been designed, constructed and 
installed, it is interesting to see, whether the wind- 
mill performs accordingly the design parameters like 
tipspeed ratio, design wind speed and efficiency and if 
not, what changes have to be made to improve the per- 
formance. The parameters which can be measured directly 
at the windmill pump system are wind speed, R.P.M. and 
output. Other parameters, like X and ^ can be derived 
from those data. 

There are in principle two types of measuring methods 
which can be applied. 

One method is measuring during a large number of very 
short periods (up to 10 seconds) during which the wind 
speed and wind direction are constant . 

These periods of constant wind speed are easy to recog- 
nize since a small variation in wind speed will result 
in an also small variation of number of revolution. 
These changes in R.P.M. are audible very well. 
The aim of this measuring method is to determine the 
tip speed - wind speed relation, efficiency - wind speed 
relation and the discharge - wind speed relation. 

Another method is measuring during longer periods: 
preferably 1 hour periods, which give more insight in 
the average performance of the windmill-pump system. 
Against hourly average wind speeds the hourly output is 
determined. This relation is more useful in output 
predictions since all losses caused by fluctuating 
wind speed and changing wind directions are included. 
Moreover, most wind data are also presented as hourly 
averages . 
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10.1 momentaneous measuring method 

The minimum equipment required for this method is: 

- windvelocity meter 

- chronometer (stopwatch) 

- small drum of about 20 1 with measuring scale 

- measuring tape to measure the elevation head 

Due to variation in wind speed, over short distance as 
well as over short periods, it is essential to have the 
windvelocity meter installed close by the rotor and in 
the same plane as the rotor. This will minimalize the 
difference of wind speed felt by the wind wheel and the 
cupanemometer. 

Accurateness is required because an error in wind speed 
reading will for instance in calculating the power 
conversion efficiency result in a deviation which is 
proportional to the cube of this error 
For measuring the amount of water pumped during the 
clocked period a 20 litre drum is held under the outlet 
of the delivery pipe during 4 or 5 strokes. The drum is 
put under the delivery during the downward movement of 
the piston when no water is being delivered. At the same 
time the chronometer is started and the strokes are 
counted. After 4 or 5 strokes the drum is shifted from 
the delivery again during the downward movement of the 
piston. At the same time the chronometer is stopped and 
the reading is printed. During the measuring cycle also 
the wind speed is printed. This complete measuring cycle has 
to be carried out with a constant wind speed. If during the 
measuring period the wind speed changes, which is very well 
detectable, the measuring cycle can be cancelled right away 
and a new cycle can be started. 

In case the wind speed increases during the measuring 
period, part of the available wind energy is used for 
accelerating the windmill-pump system, so it is not known 
how much energy is utilized for lifting the amount of 
water measured. 

If on the other hand the wind speed decreases during the 
measuring period, not only the available wind energy is 
used for lifting the water, but also part of the energy 



stored in the windmill system, which will be released by 
decelerating the windmill-pump system, so also here it 
is not known how much energy is utilized for lifting the 
amount of water measured . That is why only under station- 
ary conditions measurements should be carried out, in 
order to obtain reliable data. 

Screening of the data, i.e. selecting only the proper 

readings according certain suppositions, will become 

easier if sufficient readings have been taken. 

On the next page a measuring form can be seen. 

On the windmill site the headings and the first four 

columns are filled. At the office the other columns can 

be processed. 

R.P.M. can be found by the following formula: 

number of revolutions x 60 = [ R .p. M .] 

time interval [s] 

The water flow can be found by: 

discharge [litres] _ [l/s] 
time interval [s] 

The volumetric efficiency if)? ) is defined as: 

discharge / stroke ^ discharge 

stroke volume ^/A d s . number of revolutions 

If the discharge is given in litres the diameter and 
the stroke of the pump should be noted in dm. 

The tipspeed ratio (X) can be derived as follows 

X = 2TTn R _ 2TTR. number of revolutions 
v wind speed . time interval 

R = rotor radius [m] 
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WINDMILL PERFORMANCE MEASURING FORM 



Site: 

Dote 

Executors: 



Total Elevationhead: 
Pump stroke: 
Piston diameter: 
Windmill type: 



Serial 
Number 

- 


Wind 
speed 
reading 

(m/s) 


Time 
interval 

(sec) 


number 
of rev. 

- 


dischargi 
(liters) 


RPM 
- 


Q 

(l/s) 


Vol 
eff. 

- 


tip 

speed 
ratio 

- 


system 
eff. 

- . 


i 




















2 




















3 




















*f 




















5 




















6 








































8 




















9 




















10 




















ii 




















iZ 




















13 




















i*i 




















IS 




















16 




















i> 




















18 




















»«» 




















20 




















21 




















22 




















2S 




















21 


















25 


















26 


















2> 
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Finally the total power conversion efficiency (Cp. ^j.^) 
is calculated as follows: 

P mech „ 9w' s ' H - Q 



x 



|ov j tt R P mech 

y w -g.H.q 

Plotting the calculated data in graphs will give a clear 
picture of the relations. 

Screening the data will be easier if they are rearranged 
with an increasing windspeed sequence. Screening itself 
is a little tricky, because one tends to cancel all 
strong deviating measuring points. 

However, wrong measurements can be detected from the wind 
speed - RPH relation. In general RPM should increase with 
increasing wind speeds. If that is the general tendency of 
the measuring points, then flukes are to be imputed to acce- 
leration or decelaration effects so these measuring points 
should be omitted. 



From the graphs and the table the following design values 
can be chequed: 

- The wind speed at which the windmill-pump system works 
optimal (with its highest efficiency) should be in 
accordance with the design speed of the installation. 

- The tipspeed-ratio should be according to the design value 
at the design wind speed. 

- The volumetric efficiency of the pump should have reason- 
able values (not below 70#). 

- The input - output (v - Q) relation should be acceptable. 

On the following pages an example has been presented. 
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WIHDHTLI. PERFOHMAHCE MEASURING FORM 



Site: Total Elevationhead: m 

Date: Pump stroke: o..u m 

Executor?: Piston diameter: O.Tj m 

Windmill type: 1 PU 500 



Serial 
Number 


Windspeed 
reading 

ro/s 


Time 

interval 
(sec) 


Number 
of rev. 


Discharge 
(liters) 


HPM 


Q 

(1/sec ) 


vol. 
eff. 


tip 

speed 

ratio 


System 
effici- 
ency 


1 




P.J 


4 


13.5 


25 


1.6 


0.80 


2.99 


0.49 


m 

i 


2.7 


6.1 


4 


14.5 


39 


2.4 


0.85 


3.81 


0.59 


I 


2.9 


6.6 


4 


13.5 


36 


2.0 


0.80 


* 3.28 


0.39 


£ 


2.9 


9.4 


4 


10.5 


26 


1.1 


0.62 


2.30 


0.22 


5* 


3.1 


6.0 


4 


13.5 


40 


2.3 


0.80 


3.38 


. 0.37 


6 


3.2 


7.5 


4 


12.0 


32 


1.6 


0.71 


2.62 


0.23 


7 


3.2 


7.9 


4 


12.0 


30 


1.5 


0.71 


2.49 


0.22 


8 


3.6 


r.5 


4 


12.5 


37 


1.9 


0.74 


2.69 


0.20 


9 


3.6 


7.2 


4 


14.0 


33 


1.9 


0.83 


2.42 


0.20 


10 


4.0 


7.6 


4 


11.5 


32 


1.5 


0.68 


2.07 


0.11 


11 


4.0 


6.0 


4 


14.5 


40 


2.4 


0.85 


2.62 


0.18 


12* 


4.3 


4.5 


4 


14.5 


53 


3.2 


0.85 


3.25 


0.19" 


13 


4.3 


6.2 


4 


13.5 


39 


2.2 


0.80 


2.36 


0.13 


14 


4.6 


5.e 


4 


13.0 


41 


2.2 


0.77 


2.35 


0.11 


15 


5.2 


5.5 


4 


14.5 


44 


2.6 


0.85 


2.20 


0.10 


16 


5.2 


5.5 


4 


13.5 


44 


2.5 


0.80 


2.20 


0.09 


"7 


5.2 


4.8 


4 


14.0 


50 


2.9 


0.83 
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10.2 Contineous measuring method 



For this method it is very convenient to have a 2 - channel 
measuring recorder which registrates hourly average wind 
speeds and hourly water discharges simultaneously. 
Prom this recordings the relation between hourly average 
wind velocity and output can be derived. 

With this measuring method it is also easy to compare the 
generated power (in terms of pumped quantities of water 
over a uniform head) per square metre swept area as a 
function of the wind speed of various windmills. 
The power per m 2 is calculated with: 



P = 




[watt/m 2 



(10.1) 
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TOOL 12 PU 500 wafer pumpirg windmill 



Rotor 



diameter: 5 m 

number of blades: 12 

design of blades: 10% arched steel sheet 

tip speed ratio: 2 

maximum power coefficient: .38 



Transmission 



crank - connecting rod - crosshead system 
crank radii: .06 f .09 and .12 m. 



Security system 



automatic system. Rotor is turned out of the 
wind by indication of auxiliary vane 



Cut-in wind speed 
depending on wind regiem 



Pump 



various diameters 3", 4", 5", 6" and 8" 
depending on elevation head and wind regiem 
Single acting reciprocating piston pump with 
centrally positioned internal airchambers 



Tower 



4 leg lattice tower with standard height of 7 m. 



Materials 



2 sizes of angle iron, 1 size of flat, common 
bolts, nuts, washers etc. ballbearings, nylon or 
bronze bushes, 1mm sheetmetal and wood 



Total weight 



Output 



approx. 450 kg 
see graphs 



Price indication $ 1000 (Indian conditions) 



Status - first prototype has been built in 1977 

- further action research in 1978 with a 
dozen of these windmills in India 

- since 1980 hundreds of these windmills 
have been built 

- drawings, building manual and technical 
report available. 



